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I. INTRODUCTION 

The study of hadron production in high energy col- 
lisions is the only way to understanding structure of 
hadrons and the dynamics of the collision process. Be- 
cause of asymptotic freedom at short distance and quark 
confinement at long distance of the strong interac- 
tion, hadron formation is a soft process which cannot 
be calculated perturbatively from the first-principle for 
the strong interactions-the quantum chromodynamics 
(QCD). Therefore some model considerations are needed 
whenever hadron formation is involved in a process. One 
of the most frequently used models for hadron produc- 
tion is the parton fragmentation model. This model was 
first proposed in studying the hadron production in e + e~ 
annihilation process In such a process, e + e~ anni- 
hilate into a virtual photon which then materializes into 
a qq pair. The produced q and q move away from each 
other and each fragments or cascades into a collection 
of hadrons (or jet). Fragmentation functions have been 
used to parameterize the non-perturbative aspects of jets 
from initiating quarks and gluons in last a few decades. 
As stated clearly in : 2f, however, this approach "is not 
meant to be a theory. It is simply a parameterization 
that incorporates many of the expected features of frag- 
mentation." The fragmentation functions does not tell 
us how the hard partons fragment into jets. If hard par- 
ton fragmentation were the only way for the production 
of hadrons with large transverse momentum, the ratio 
between two hadron species, for any given transverse 
momentum, would be a process independent constant 
dictated by the fragmentation functions. However, ex- 
perimental data at the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven have shown that the proton to 
pion ratio can be about 1 for pt ~ 3 GeV/c in central 
AuAu collisions at ^/snn — 200 GeV Q, much larger 
than the corresponding value in pp collisions. This exper- 
imental fact indicates that the study of the mechanism of 
parton fragmentation is necessary in order to understand 
the colliding system dependence of hadron production in 
a jet. 

A theoretical microscopic description of parton frag- 
mentation is suggested in Q based on quark recombina- 
tion approach |5|, |6j for hadron production. The physics 
issue in the microscopic description of fragmentation is 
that every quark and gluon with high virtuality is as- 



sumed to evolve by emanating quark-antiquark pairs and 
gluons with lower virtualities into a shower of partons 
with low transverse momentum relative to the initiating 
parton direction. The shower partons recombine then to 
form the final state hadrons. In this way, the fragmen- 
tation functions can be related to the shower parton dis- 
tributions in a jet. In high energy hadronic and nuclear 
collisions, there exist some soft partons with low virtu- 
ality and some hard partons produced in hard collisions 
with high momentum transfer. Those hard partons may 
have very high virtuality and, according to the principles 
of QCD, will evolve into parton showers. Those shower 
partons, together with the soft partons, then recombine 
into final state hadrons detected in experiments, as de- 
scribed in the recombination model. 

In Q shower parton distributions initiated by a hard 
parton with virtuality Q 2 — 100 (GeV/c) 2 were obtained 
from the fragmentation functions determined by fitting 
experimental data for high energy e + e~ annihilation and 
hadronic collisions. Once the shower parton distributions 
are obtained, it is natural to apply them in modeling 
hadron production in other collisions. In 6] it is shown 
that the recombination of soft partons and shower par- 
tons is important in reproducing the transverse momen- 
tum spectra of produced particles in both AuAu and dAu 
collisions and in explaining the centrality dependence of 
the Cronin effect pj in dAu collisions. 

In this paper, we investigate the parton distributions 
in a jet for Q 2 from 2 to 1000 (GeV/c) 2 . We also param- 
eterize the Q 2 dependence of the parameters for the par- 
ton distributions for easy use in the recombination model 
description for hadron productions in different collision 
processes. The application of the parton distributions 
will be left for later work. This paper is organized as 
follows. In section II, we briefly discuss the recombina- 
tion model description of parton fragmentation. Then 
in section III, we discuss the Q 2 evolution of the parton 
distributions. The last section is for a short discussion. 



II. SHOWER DISTRIBUTIONS IN A JET 
FROM RECOMBINATION MODEL 

For brevity, one can introduce following notation for 
any two parton distributions fi(x) and f2(x) in a jet, 

{fif2}(xi,x 2 )= 0.5(/i(a:i)/a( T ak r ) + 
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Here xi,X2 are momentum fractions of two partons rel- 
ative to the initiating parton. In this expression, kinetic 
constraint is imposed so that x\ + X2 < 1 for any pair of 
two quarks. 

In the quark recombination model for fragmentation 
the generic formula for the production of a meson com- 
posed of a quark q and an antiquark q' is 0, @ 

xD{x) = [ — / < ^{qq'}{xx,x-z)R(xi,X2,x) , (1) 

Jo x l Jo x 2 

where q{x\) and (f (2^2) are the parton distributions of 
a quark q at momentum fraction x\ and an antiquark 
q' at X2 in a jet, and R(x±, X2, x) is the recombination 
function (RF) for the formation of a meson at x from the 
quark q and antiquark q' . In the recombination model, 
the gluons are assumed to convert first to quarks and 
antiquarks before hadron formation. So in last formula 
no gluon distribution is considered. 
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FIG. 1: Fitted fragmentation functions for valence quark to 
pion, sea quark to pion, gluon to pion, and sea quark to kaon, 
gluon to kaon, at (a) Q 2 = 10 and (b) 100 (GeV/c) 2 . The 
fitted results are shown as solid curves, and the parameteriza- 
tions from KKP (for gluon jets) and Kretzer (for quark jets) 
by different symbols for different fragmentation functions. 

At this stage, we only consider production of hadrons 
composed of light quarks (m, u, d, d, s, s) fragmented from 
light hard partons (u, u, d, d, s, s, g). If we assume SU(2) 
symmetry and symmetric uu and dd sea quark distri- 
butions, there are five independent parton distributions. 
Three of them are for parton distributions in a jet ini- 
tiated by light quarks. They are the non-singlet distri- 
bution i^NS of valence quark with the same flavor as the 
initiating (anti)quark, the non-strange sea quark distri- 
bution L, the strange sea quark distribution L s . The 



other two are for the non-strange light quark distribu- 
tion G and strange quark distribution G s in a gluon jet. 
Five fragmentation functions are needed in order to de- 
termine those five parton distributions. We work with 
fragmentation functions for the pion production D§ from 
sea quarks, D v from valence quark, D G from gluon jet, 
and for K production Dg from sea quarks, Dq from 
gluon jet. The corresponding expressions for the frag- 
mentation functions in terms of the parton distributions 
are given in the quark recombination model as Q 



xDl{x) 



xD v {x) 



xD G (x) 



xD§(x) 



dxi dx2 



Xl X2 



{LL}(x 1 ,x 2 )R w (x 1 ,x 2 ,x), (2) 



dx\ dx2 



Xl x 2 



dxi dx2 
Xl x 2 



{KL}(xi,x 2 )R 7T (xi,X2,x), (3) 



{GG}(xi,x 2 )R v {x u x 2) x), (4) 



dxi dx 2 
Xl x 2 



{LL s }(xi,x 2 )R K (xi,X2,x), (5) 



xD%(x)= j ^^-{GG s }(xi,x 2 )R K (xi,x 2 ,x). (6) 

J Xl X 2 

In last equations i? T and R K are RF for pions and kaons, 
respectively. Here, both the fragmentation functions and 
the parton distributions in a jet depend on the virtuality 
Q 2 of the initiating parton. 

Before one can start determining the parton distribu- 
tions from above equations, one needs to make a decision 
on choosing parameterizations of the fragmentation func- 
tions given by different groups. In fact, the determina- 
tion of fragmentation functions is not as well as supposed 
to be. Different physics considerations are input in dif- 
ferent parameterizations of the fragmentation functions, 
and sometimes the obtained results are not consistent 
with each other. The most frequently used parameteriza- 
tions of the fragmentation functions are given by Kniehl, 
Kramer, Potter (KKP) [1G3 and by Kretzer (Tj. In [HI 
it is assumed that valence quark dominates the produc- 
tion of hadrons with large momentum fraction and sea 
quarks contribute more to hadron production at small 
momentum fraction. In KKP parameterizations, no such 
physics consideration is input, but they used more re- 
cent experimental results involving gluons. In our paper, 
we adopt Kretzer parameterization of the fragmentation 
functions for jets initiated by quarks or antiquarks for a 
better consideration of valence quark contributions, but 
use KKP parameterization for gluon jets. 

In fitting the fragmentation functions from quark re- 
combination model, as given in Eqs. ill'Kil) . we assume 
the parton distributions to be the form 



f(z) = Az a (l-z) b (l + cz d ) , 



(7) 
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with five parameters A, a, b, c and d, and z the momen- 
tum fraction of the quark relative to the initiating par- 
ton. For the valence quark distribution, a condition is 
imposed 



dz 



(8) 



because the total number of valence quark in a jet cannot 
be more than one. 

The fragmentation functions for pion and kaon have 
been fixed in 113 as 



n*t \ X1X2 [ X1+X2 . 
R (Xl,X2,X) = — 5-d I 1 



(9) 



R K (xi,x 2 ,x) = 12^5 (£i±£a _ A . (10) 

X° \ X J 

In the recombination function R K for kaons, x 2 refers to 
the momentum fraction of the strange (anti)quark and 
x\ to that of the non-strange (anti)quark. 

Then one can fit the fragmentation functions and get 
the parton distributions in a jet for any fixed Q 2 . The 
obtained parameters A, a, b, c and d for the five parton 
distributions depend on the scale Q 2 . For Q 2 = 10 and 
100 (GeV/c) 2 the fitted results are shown in Fig. 1 as 
solid curves, in comparison with the KKP and Kretzer 
parameterizations. Obviously, the recombination model 
can fit the fragmentation functions very well. It should 
noted that the fitted curve for K-^s agrees with Kret- 
zer parameterization quite well, in sharp difference from 
the fitted result in |4| , there the parameterization for the 
valence quark fragmentation had a significant contribu- 
tion at small momentum fraction region, while in Kretzer 
parameterization the valence quark fragmentation is as- 
sumed to dominate in large momentum fraction region 
and no singularity at very small momentum fraction. 

The fitted shower parton distribution functions 
(pdf(z))&re shown in Fig. 2 for Q 2 = 10 and 100 
(GeV/c) 2 . As expected, the non-singlet distribution K^s 
for valence quark is hardest among the five distribu- 
tions. Other parton distributions show a strong increase 
at small z as z goes to smaller values. 



Q. 




FIG. 2: Fitted shower parton distribution functions (pdf(z)) 
in a jet at (a) Q 2 = 10 and (b) 100 (GeV/c) 2 . 
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III. EVOLUTION OF PARTON 
DISTRIBUTIONS 

In only the fragmentation functions at fixed Q 2 = 
100 (GeV/c) 2 are employed. The obtained parton dis- 
tributions are then used to calculate the contributions 
from both soft-shower and shower-shower parton recom- 
bination to the hadron px spectra. For simplicity, all 
hard partons are assumed to have the same virtuality 
Q 2 = 100 (GeV/c) 2 in all applications. In more realistic 
cases, partons with different momenta may have differ- 
ent virtualities. This difference in virtuality may influ- 
ence hadron productions from the fragmentation of hard 
partons and from other quark recombination processes. 



FIG. 3: Q 2 dependence of five parameters for parton distribu- 
tions -Z^ns and L. Solid curves are fitted results by Eq. 111^ . 
The QCD scale A is chosen to be 0.19 (GeV/c) 2 . 



From perturbative QCD we know that the Q 2 evolu- 
tion of the fragmentation functions is described by the 
DGLAP equations 13]. From such equations one can 
determine the fragmentation functions at any Q 2 if one 
knows them at one given scale Qq. Up to now however, 
we did not know the Q 2 dependence of parton distribu- 
tions in a jet initiated by a hard parton, nor have we had 
an equation to describe their evolution in a model inde- 
pendent way. But we can get such a dependence by using 
Eqs. (I-?l()l) from the Q 2 evolution of the fragmentation 
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functions. In this paper, we simply use the FORTRAN 
code provided by Kretzer 01 to get fragmentation func- 
tions at any given Q 2 . For each Q 2 we use Eqs. (12101) to 
fit the fragmentation functions and get the correspond- 
ing parameters for each parton distributions at that scale 
Q 2 . 

In accordance with other theoretical results, we define 
here a dimensionless scale u = ln(Q 2 /A 2 ) with the QCD 
scale A = 0.19GeV/c. It is found that the parameters for 
each shower parton distribution has a strong dependence 
on the scale u. As examples, the Q 2 dependence of pa- 
rameters A, a, b, c, d for L and i^NSi obtained from fitting 
the fragmentation functions at different Q 2 , are shown in 
Fig. 3 as functions of u = \n(Q 2 /A 2 ) by symbols. For 
easy use of others, one can express the Q 2 dependence of 
the parameters by a polynomial of the scale u as 

A(Q 2 ) = C u 3 + C lU 2 + C 2 u + C 3 . (11) 

The fitted results are shown by solid curves also in Fig. 3 
for comparison. Then the dependence of every of five pa- 
rameters for a shower parton distribution in a jet is given 
by four parameters, Co, Ci, Ci and C3 in a wide range of 
Q 2 from 2 to 1000 (GeV/c) 2 . Since the Q 2 dependence 
of each parameter in the parton distributions can be well 
reproduced by only four constant coefficients, Eq. I|llfl is 
a quite economic description of the evolution of shower 
parton distributions. The coefficients Cq,Cx,C% and C3 
are tabulated in Table I for all the parameters for the 
five parton shower distributions discussed in this paper. 



From this table, parton distributions can be calculated 
easily for a jet with any scale Q in the range considered 
in this paper. 



IV. DISCUSSIONS 

The evolution of shower parton distributions in a jet is 
discussed in the framework of quark recombination model 
from the Q 2 dependence of the fragmentation functions. 
We show that the evolution of the distributions can be 
given by the polynomial InQ 2 dependence of a few pa- 
rameters. From the tabulated coefficients one can cal- 
culate corresponding parton distributions in a jet and 
other fragmentation functions needed in the framework 
of quark recombination model. Together with soft and 
hard parton distributions, these shower parton distribu- 
tions can be used for the production of hadrons in differ- 
ent collision processes. 
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TABLE I: Fitted coefficients in Eq. 11 U for parameters of 
the parton distributions in a jet. 







Co 


Ci 


c 2 


c 3 


Kns 


A 


2.27 x 1(T 5 


-4.91 x 10~ 4 


2.78 x 10~ 3 


5.92 x 10" 3 




a 


-7.21 x 1(T 5 


1.94 x 10~ 3 


-1.92 x 10~ 2 


8.46 x 10~ 2 




b 


1.24 x 1(T 3 


-3.47 x 10" 2 


0.398 


-6.24 x 10~ 3 




c 


-0.328 


7.67 


-45.3 


246 




d 


2.92 x 1(T 4 


-7.17 x 10~ 3 


4.88 x 10~ 2 


0.993 


L 


A 


-0.0009 


0.0275 


-0.3179 


1.7997 




a 


-0.0003 


0.0116 


-0.1489 


0.4376 




b 


0.0023 


-0.0588 


0.5344 


1.6756 




c 


0.0251 


-0.5642 


3.8456 


-2.7196 




d 


0.0040 


-0.0773 


0.2938 


5.2021 


L a 


A 


-0.0427 


0.9992 


-7.7046 


21.3488 




a 


-0.0049 


0.1210 


-1.0107 


3.5884 




b 


-0.0022 


0.0442 


-0.1268 


8.1004 




c 


0.1677 


-4.2572 


38.9198 


-18.9039 




d 


-0.0014 


0.0308 


-0.2036 


3.4815 


G 


A 


-0.0025 


0.0835 


-0.9312 


3.9350 




a 


0.0004 


0.0007 


-0.1318 


0.5742 




b 


-0.0015 


-0.0030 


0.4695 


-0.2050 




c 


0.0003 


-0.0073 


0.0575 


-1.1224 




d 


0.0014 


-0.0906 


1.0256 


-1.6076 


G s 


A 


-0.0015 


0.0319 


-0.2183 


0.5004 




a 


-0.0140 


0.3064 


-2.0830 


3.6414 




b 


-0.0061 


0.1095 


-0.3339 


1.6614 




c 


0.1277 


-2.7241 


18.0042 


-34.0906 




d 


-0.0128 


0.2812 


-1.9713 


5.6142 



